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ABSTRACT 

We present high-sensitivity multi-frequency radio continuum observations of the star- 
burst galaxy NGC 2146. We have fitted these data with a three-dimensional diffusion 
model. The model can describe the radio emission from the inner disk of NGC 2146 
well, indicating that diffusion is the dominant mode of propagation in this region. Our 
results are indicating that NGC 2146 has recently undergone a starburst, the star form- 
ing activity being located in a central bar. The spatial variation of the radio emission 
and of the spectral index yield tight constraints on the diffusion coefficient Do and the 
energy dependence of the diffusion. Away from the central bar of the galaxy the radio 
emission becomes filamentary and the diffusion model was found to be a poor fit to the 
data in these regions; we attribute this to different transport processes being important 
in the halo of the galaxy. 
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1 INTRODUCTION 

Multifrequency radio observations of nearby, face-on galax- 
ies present a unique way to study the Cosmic Ray (CR) 
propagation within their disks. These studies provide con- 
straints on the main parameters of this propagation, i.e. the 
value and energy dependence of the diffusion coefficient and 
the form of the electron energy distribution immediately af- 
ter acceleration. Furthermore, because of the relatively long 
radiative life-time of CR electrons (which is in normal, non- 
starburst galaxies up to several 10 7 years) such observations 
also provide information about the temporal evolution of the 
sources of the CR's on the same time scale; in turn, since 
the sites of acceleration are supernova remnants the data 
are probing the history of star formation on time scales up 
to of order « 10 8 years. 

To date, the parameters describing CR propagation 
have only been derived for our Galaxy; radio observations 
are the only way of obtaining information for external galax- 
ies. In this paper we analyse high-sensitivity multi-frequency 
radio data for the nearby galaxy NGC 2146. Our data have 
sufficient frequency coverage that we can include in our 
model fits both a thermal and diffuse synchrotron compo- 
nent for the radio emission. 

NGC 2146 is a very luminous radio 
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GHz = 5.910^ W/Hz) at a distance of 21.; 
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Mpc (Benvenuti, Capacioli & D'Odorico 1975, adjusted to 
-ffo=50 km s _1 Mpc). Most of its radio luminosity is due to 
a compact starburst in its centre. The extent of the radio 
emission (1.5 arcmin « 9.6 kpc) is much less than the extent 
of the optical emission (« 6 arcmin) . This suggests that ac- 
tive star formation (both current and in the recent past) is 
confined to this central region of NGC 2146. In this paper 
we are only discussing this inner region of the galaxy (w 1.5 
arcmin x 0.7 arcmin) that is characterized by the strong 
radio emission. Optical images of NGC 2146 show it to be a 
peculiar spiral galaxy (e.g. Benvenuti et al. 1975). Hutchings 
et al. (1990) have presented a detailed study of NGC 2146 
in several wavebands and proposed that the galaxy is in the 
late stages of a merger which has given rise to an intense 
burst of star formation in the nucleus. 



2 OBSERVATIONS 

NGC 2146 was observed with the VLA at 8.4 GHz and 
327 MHz, and with the Cambridge 5-km Ryle Telescope 
(RT) at 5 and 15 GHz. Additionally, existing VLA data 
at 1.5 GHz (Condon 1983) were available. When the ob- 
servations of NGC 2146 were made, the RT was operated 
with a band-width of 280 MHz split into 28 10-MHz fre- 
quency channels; together with the minimum baseline of 
18 m this provides excellent temperature sensitivity at both 
5 and 15 GHz. 

For the RT observations a total of five telescope config- 
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Figure 1. Maps of NGC 2146 at 1.5, 5, 8.4 and 15 GHz at a resolution of 4.7 X 4.7 arcsec 2 . The contour levels 
start at 1 mjy/beam (1.5GHz) and 0.5 mjy/beam (for the other frequencies), and increase by factors of \f 7 l. The 
r.m.s. noise levels are: 0.67 mjy/beam (1.5 GHz), 0.244 mjy/beam (5 GHz), 0.18 mjy/beam (8.4 GHz) and 0.16 
mjy/beam (15 GHz). 



urations at each observing frequency were employed giving 
a nearly fully filled aperture out to the maximum baseline. 
A minimum of two 12-hour runs in each configuration were 
obtained; any data of poor quality were rejected. This gave a 
resolution approximately 1.0 x 1.0cosec(<5) arcsec 2 at 15 GHz 
and 3.0 x 3.0cosec(<5) arcsec 2 at 5 GHz. At 5 GHz observa- 
tions of phase-calibrators were made at the beginning and 
end of each run, while at 15 GHz calibration observations 
were interleaved with those of NGC 2146. 3C 286 and 3C 48 
were observed regularly as flux calibrators. Calibration and 
data-editing (principally to remove narrow-band and time- 
varying interference) was performed in the MRAO package 
Postmortem, with subsequent reduction in AIPS and the 
MRAO package Anmap. 

The 8.4-GHz VLA observations were in the B, C and 
D arrays with the 327 MHz observations in A array only; 
the observations were performed during 1990/91. Reduction 
of the data followed standard VLA procedures with calibra- 
tion, editing and imaging performed in AIPS. 

A number of maps were made at each frequency by ta- 
pering the aperture-plane to highlight structure on different 
angular scales. For the spectral comparison we used only 



the four highest frequencies and all the maps were made 
with very similar aperture plane coverage, CLEANed and 
restored (or convolved in the case of the 1.49 GHz image) 
to a common resolution of 4.7 x 4.7 arcsec 2 . All maps have 
been corrected for the primary beam response of the tele- 
scope used to make the observations. This set of four images 
is shown in Fig. 1, together with a high-resolution image at 
15 GHz in Fig. 2 and the 327 MHz image in Fig. 3. In Fig. 
4 a map of the spectral index between 1.5 and 8.4 GHz is 
shown. 



3 RADIO PROPERTIES OF NGC 2146 

Most of the radio emission of NGC 2146 lies along an elon- 
gated, bar-like structure in position angle 143 ± 3° coin- 
cident with the optical major axis (143°, Benvenuti et al. 
1975). This elongated shape could be a disk viewed from a 
high inclination or it could be due to a real bar-like shape 
of the sources of the CR's. From the radio data at one fre- 
quency alone it is not possible to distinguish between these 
two cases. However, as will be shown in the following, the 
variations of the radio spectral index together with a model 
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Figure 3. Map of NGC 2146 at 327 MHz at a resolution of 8 
X 8 arcsec 2 . The contour levels are start at 10 mjy/beam, and 
increase by factors of y/2. The r.m.s. noise level is at 8 mjy/beam. 



for the CR propagtion allow to draw conclusions about the 
shape of the CR sources. In the high resolution image (Fig. 
2), the sources at 15 GHz show a distinctive S shape that 
had already been noticed in earlier observations by Kron- 
berg & Biermann (1981). The central region is shown in our 
highest resolution images to contain a number of luminous 
unresolved sources superimposed on the general background 
emission. 

At higher frequencies NGC 2146 appears more compact 
and shows relatively little emission away from the major axis 
of the radio emission while at lower frequencies the emission 
has an irregular appearance and several extensions perpen- 
dicular to the major axis which are prominent particularly 
in the 1.5 GHz and the 327 MHz images. 



4 MODELLING THE RADIO EMISSION 

Our aim in this work is to explain the multi-frequency radio 
data by a combination of thermal emission and non-thermal 
synchrotron radiation from CR electrons which have propa- 
gated away from the sites where they have been accelerated. 
The thermal component arises in Hn regions, whereas su- 
pernova remnants (SNRs) are believed to be, at least from 
energetic considerations, the most likely source of CRs (e.g. 
Berezinsky et al. 1990, Volk 1992). 

The propagation of the CR electrons is described by a 
standard three dimensional diffusion model (e.g. Ginzburg 
& Syrovatskii 1964). Other modes of propagation, such as 
convection, may also be important especially in the halo 
(e.g. Hummel & Dettmar 1990) at large distances from the 
galactic plane, but for the sake of simplicity we do not take 
this into account explicitly. Our approach, therefore, is to 
determine to what extent we can explain the observations 
using a simple diffusion model alone. 

The diffusion equation describing the evolution of the 
electron distribution function f(E, r, t) is as follows: 



df ^ + D(E) V \f(E,r,t) 



dt 



- — [b (E)f(E,r,t) 



-Q{E,t,v) 



(1) 



where D(E) is the diffusion coefficient, assumed to be spa- 
tially homogeneous and isotropic. The diffusion coefficient 
is assumed to depend on the energy as 



D(E) = Do 



E 



lGeV 



(2) 



which is in agreement with observations of CRs in our 
Galaxy (see e.g. Berezinsky et al. 1990). In the energy range 
relevant to our radio observations (E > 1 GeV), the energy 
losses of the electrons are mainly synchrotron and inverse 
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Compton 
b(E) = 



,dE 



_dE 
dt 
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(f/rad + f/e) 



(3) 



where <tt is the Thompson scattering cross section, c the 
speed of light, ?7b the energy density of the magnetic field 
and C/rad the energy density of the radiation field below the 
Klein-Nishima limit (e.g. Longair f992). 

For the source function Q(E, r, t) we make the following 
assumptions: 

(i) The total source distribution can be separated and 
expressed as the sum of individual sources i.e. 



(4) 



(ii) CRs are produced only by SNRs from massive stars 
(m > 5Mq). The source spectrum is assumed to be a power- 
law in energy and to be identical for all sources 



£i{E) = 



^SN<?SN- 



(5) 



Here, vsn is the SN rate and qsN is the number of relativistic 
electrons produced on average per energy interval by a SNR. 

(iii) The spatial distribution of a source is approximated 
by an ellipsoid, gaussian distribution 



fci(r) 
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with (xi, yi,Zi) indicating the site of the centre of the source 
and R x ,iRy,iRz,i the axes of the ellipsoid. In the following, 
we will identify the major axes of NGC 2146 with the x-axis, 
the direction within the galactic plane perpendicular to the 
major axis with the y- and the direction perpendicular to 
the plane with the z-axis. 

Eq.(l) can be solved for an infinite space following the 
analysis of Syrovatskii (1959). In a steady-state situation 
(i.e. for constant sources) the electron distribution function 
is given by: 

f &2 

f(E, r) = vsNqsN 



exp 
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where E2 , E\ are the upper and lower end of the range of 
electron injection energies and A is given by: 

' ,;o ^' w 
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c(i-m)"" v " "° r (8) 

If temporal changes are important, a similar solution 
for f(E, r, t) can be derived. The synchrotron emission 
Psyn{v, r, t) is derived from f(E, r, t) with the simplifying as- 
sumption that each electron emits its synchrotron radiation 
at the maximum frequency of the synchrotron spectrum. 



P svn {v,r,t) = f{E(v),r,t) — 

with the emitting frequency being 

E - 2 
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(9) 



(10) 



where vg = eB/m e c 2 is the gyrofrequency. 

The thermal emission is also described as arising from 
ellipsoidal sources, that are not necessarily the same as the 
CR sources 



Ptherm(v, r) OC 



(11) 



Both, P 3 yn and Ptherm have to be integrated through the 
halo along the line of sight, taking into account the inclina- 
tion angle of the galaxy, and convolved with the beam. The 
total radio emission, P ra d, is finally obtained as the sum of 

Psyn and Ptherm- 



5 FITTING THE MODEL TO THE RADIO 
DATA 

In order to keep the diffusion model as simple as possible we 
shall assume a steady state, i.e. we attempt to fit a model 
in which the star formation rate and hence the sources of 
thermal and non-thermal emission do not show any tempo- 
ral variations. In NGC 2146 the estimated life-time of CR 
electrons is, due to the high radiation losses (see below), ex- 
tremely short (less than 2 x 10 6 years). On these short time 
scales the assumption of a steady state seems reasonable. 
Furthermore we will assume as the simplest case that Hn 
regions and supernova remnants are co-spatial (at least at 
the resolution of our radio images) . This will be a good ap- 
proximation if the sites of star formation have not changed 
during the past ~ 10 8 years, i.e. the life-time of supernova 
progenitors (m > 5M ) . If there has been significant propa- 
gation of the star formation process through the galaxy then 
this assumption must be modified. 

Additionally we can constrain some of the parameters 
directly from observational data. 

1) The inclination of NCG 2146 is = 65° (Benvenuti 
et al. 1975). Due to the disturbed appearance of NGC 2146 
the inclination angle is unlikely to be precisely determined; 
we shall adopt this value and discuss this point more fully 
in Section 6.1. The radius of the CR sources perpendicular 
to the disk is taken to be R Zii — lOOpc, which is a typ- 
ical scale height of the molecular gas disk (e.g. Unsold & 
Baschek 1991) in which we expect massive star formation to 
be located. 

2) We estimate the magnetic field B by applying a mini- 
mum energy calculation. We adopt standard values for most 
parameters: a frequency range of 10 7 -10 10 GHz; radio spec- 
tral index 0.7 (here and in the following we define the spec- 
tral index a as P ra d <x v ~ a )\ ratio of total energy of CR's to 
total energy of electrons k = 100. For measured quantities 
we use the emission at 1.5 GHz. We obtain the volume of 
the emitting region at this frequency by assuming that it 
can be approximated by a cylinder with a height 1 kpc and 
a radius of 45 arcsec (= 5kpc) which corresponds to the low- 
est contour of the 1.5 GHz image. With these assumptions 
we obtain B — 38^iG. 
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Figure 5. Slices of the radio emission at 1.5, 5 and 8.4 and 15 GHz along the minor and the major axis (full 
lines), together with the best-fit steady state model results (dashed lines). The parameters are: 7 = 1.7, fi = 0, Do = 
2 X 10 28 cm 2 s _1 , P therm (1.5GHz) = 4 mjy/beam at r = 0, R x a = 1.3 kpc, R yA = 0.7 kpc, R Xt2 = 3.0 kpc, 
Ry,2 = 0.5 kpc. 



The energy density of the radiation field, [Trad, can be 
calculated from the bolometric luminosity: 



U„ 



2L bol _ (Ltoi/W) 
(i?/kpc) 2 



cnR 2 



1.4 • 10" 



eV 



(12) 



For NGC 2146 the greatest part of its luminosity is emitted 
in the far-infrared (FIR); we therefore approximate the bolo- 
metric luminosity L bo i by the far- infrared luminosity Lfir- 
For the radius R we take the radius derived from our ra- 
dio data (R — 5 kpc) as this will correspond to the region 
in which most star formation is currently occurring. With 
these assumptions we obtain U ia ,d = 20 eVcm~ 3 . The en- 
ergy density of the radiation field, U ra d, can be calculated 
from the bolometric luminosity: 



U„ 



2L bol _ (Lboi/W) 
(R/kpc) 2 



cnR 2 



1.4 • 10" 



eV 



(13) 



For NGC 2146 the greatest part of its luminosity is emitted 
in the far-infrared (FIR); we therefore approximate the bolo- 
metric luminosity L bo i by the far- infrared luminosity Lfir- 
For the radius R we take the radius derived from our ra- 



dio data (R — 5 kpc) as this will correspond to the region 
in which most star formation is currently occurring. With 
these assumptions we obtain f/ ra d = 20 eVcm -3 . 

With these derived values, the fitting of the diffusion 
model now proceeds by fitting to slices extracted from our 
mult i- frequency radio images along both the major and mi- 
nor axis of NGC 2146. In Figs. 5 and 6 we show the obser- 
vational data and the corresponding model fits. 

In constraining the parameters for the model we note 
some general features of the data. 

i) The spectral index between 1.5 and 8.4 GHz is ap- 
proximately constant along the major axis (Fig. 8a); this 
suggests that the CR source is itself extended in this di- 
rection unless the diffusion coefficient had a strong energy 
dependence with /x « 1. (In this latter case A (Eq. 8) would 
be largely independent of the electron energy and hence the 
spatial distribution of the electrons would also be energy 
independent (see Eq. 7)). 

ii) The steepening of the spectrum away from the major 
axis (Fig. 8) indicates that fj, is considerably less than one 
and that the electrons must have suffered from energy losses. 
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iii) The radio flux decreases away from the centre of the 
galaxy both along the major and minor axis (Figs. 5 and 6). 

Taken together, these observations suggest that the 
model should have the following general characteristics: 

• The source distribution is extended along the major 
axis but not along the minor axis, i.e. it has the shape of a 
bar. 

• The energy dependence of the diffusion coefficient, /i, 
must be considerably less than unity, and is constrained by 
the data along the minor axis of NGC 2146. 

The fitting of the model to the observational data now 
proceeds by a straightforward search of parameter space for 
model parameters which lead to the best fit of the multi- 
frequency radio data. The model most strongly constrains 
the spatial distribution and strength of the thermal and non- 
thermal emission regions, and also the energy dependence of 
the diffusion coefficient. Our results and conclusions for the 
simplest form of the model without temporal variations are 
shown in Fig. 5 and can be summarized as follows. 

1) The source structure cannot be adequately repre- 
sented by a single ellipsoidal emission region. Two sources, 
both centred at the maximum of the radio emission, are nec- 
essary to account for the data. Our best fit parameters for 
these two sources are: 

- source 1: R Xt i — 1.3 kpc, R y> i — 0.7 kpc, relative inten- 
sity 1 

- source 2: R Xt 2 = 3.0 kpc, R y> 2 = 0.5 kpc, relative inten- 
sity 0.3 

We interpret this result as suggesting that the distribution 
of CR sources and hence star formation activity is in the 
form of a bar. The parameters describing source 1 are rel- 
atively well constrained by the shape of the radio emission 
(for a fixed value of the galaxy inclination). The parame- 
ters of source 2 are less well constrained, however this has a 
relatively minor effect on the other model parameters. 

2) The injection index of the electron population is 7 = 
1.7 and we infer a low thermal fraction of Ptherm/ Prad = 17 
per cent at 15 GHz. The thermal fraction is constrained by 
the radio spectrum at the centre of the galaxy. 

3) It is clear from Fig. 5 that we are not able, with the 
simplest form of the model, to get good fits to the data over 
the whole frequency range: At 5, 8.4 and 15 GHz the model 
predictions for the peak radio emission are either slightly too 
low or too high. The discrepancies however never exceed 10 
% and it should furthermore be noted that the shape of the 
radio emission can be reproduced well at all frequencies. 

The results obtained by allowing for temporal varia- 
tions in the star formation rate, while keeping the source 
distribution constant are shown in Fig. 6. In this case, an 
injection index of 7 = 1.65 is used and the temporal vari- 
ations are assumed to have a simple form with injection of 
CR electrons having stopped 5.7 x 10 5 years ago. The pre- 
dicted thermal contribution is now much greater than in the 
steady-state model. At 15 GHz the emission is almost com- 
pletely thermal and at 1.5 GHz the thermal fraction is still 
Ptherm/ P ra d = 27 per cent. The time-dependent solution 
produces an improved model fit as it permits an aged elec- 
tron population to develop. In general, for both the steady- 
state and time-dependent cases, the fits along the major axis 
are worse than those along the minor axis — this is because 



of the structure of P ra d which reflects the complicated in- 
trinsic source structure (see Fig. 2) which we have not taken 
into account. However, this does not significantly effect are 
conclusions concerning propagation of CR to be discussed 
below. 

4) The value of the diffusion coefficient Do , and its en- 
ergy dependence, is strongly constrained by the data per- 
pendicular to the inferred bar within NGC 2146 along the 
minor axis. Furthermore, the results are very similar for both 
the steady-state and a time-dependent models. In Fig. 7 we 
show the results for the steady-state case - the best fit for 
the diffusion coefficient is D = 1.0 -> 4.0 x 10 28 cm 2 s _1 , 
with [i, = 0. The shape of the spectral index distribution 
along the minor axis depends sensitively on /x, which is il- 
lustrated in Fig. 8. Our best estimate for fi is in the range 
/j, = 0.0 — » 0.2 with a firm upper limit of pi < 0.5. 



6 DISCUSSION 

6.1 Model uncertainties 

As illustrated in Fig. 5-8 the model describes the data well 
and allows us to constrain important parameters such as the 
diffusion coefficient and its energy dependence and infer the 
existence of a central bar. The error in the determination of 
these parameters depends however not only on the quality of 
the model fit, but also on a number of parameters that were 
determined independently such as B,U T!l d, the inclination 
angle and distance of the NGC 2146. 

In order to fit the spectral data a model that produces 
a break in the synchrotron spectrum is required - this is 
evident from the data since the spectrum at the centre of 
NGC 2146 steepens between 5 and 8.4 GHz before flattening 
towards 15 GHz which we attribute to thermal emission. 
The simplest way to achieve such a break is to assume some 
temporal dependence of the CR injection and hence of the 
star- formation rate. The time-scale of the changes implied by 
the data is rather short, less than 10 6 years, however there is 
a strong dependence on B and U r3 A'- t oc B 1 ^ 2 / {U IS ,d + Us). 
Whereas our estimate for (7 rat j does not suffer from large 
uncertainties, B is very uncertain and an improved estimate 
would require an independent determination of the magnetic 
field. 

The main sources of error for the diffusion coefficient, 
Do, are the uncertainty in [7 ra d, B, the distance to NGC 2146 
and the assumed structure of the emitting region. In princi- 
ple it could be possible to increase/decrease the size of the 
CR source along the minor axis (R y ) and try to fit the data 
by decreasing/increasing Do- If R y is very small, even with 
a large value of Do the model cannot produce the observed 
shape of P ra d- On the other hand, a larger assumed source 
size would make it impossible to fit the spectral variations 
unless we assumed that either the injection index of the CR 
population varied within the disc of the galaxy, and/or that 
different regions of the galaxy had undergone very differ- 
ent star-formation histories. Within the range of acceptable 
R y the estimated values for Do all lie within our quoted 
error bounds. The uncertainty in the inclination angle of 
NGC 2146 does not affect our estimates of Do and \x greatly 
and it most strongly affects the assumed source size. Our 
conclusion regarding the bar-like shape of the star forming 
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Figure 6. Slices of the radio emission (full lines), together with the best-fit results for a model with temporal changes 
in the injection rate (dashed lines). The parameters are: 7 = 1.65, fj, = 0, Do = 2 X 10 2s cm 2 s _1 , P therm (1.5GHz) = 
23 mjy/beam at r = 0, the source structure is the same as in Fig. 5. The injection of CR's has stopped 5.7 X 10 s 
years ago. 




Figure 7. Slices of the radio emission at 1.5 GHz along the minor and the major axis (full line) together with model 
results for different values of the diffusion coefficient: Do = 1.0 X 10 28 cm 2 s~ 1 (dashed line), Do = 2.0 X 10 2s cm 2 s~ 1 
(dashed-dotted line), and Do = 4.0 X 10 28 cm 2 s _1 (dotted line). The other parameters are as in Fig. 5. 
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Figure 8. Slices along the minor and the major axis of the spectral index between 1.5 and 8.4 GHz, ct(i.5_8.4G-Hz) 
(full line), together with the model results for different values of fi. The parameters are: \x = 0,7 = 1.7, Do = 
2.0 X 10 28 cm 2 s _1 (dashed line), \i = 0.2,7 = 1.65, A) = 1.5 X 10 28 cm 2 s _1 (dashed-dotted line), fi = 0.5,7 = 
1.6, D = 0.8 X 10 28 cm 2 s _1 (dashed-dotted line), and fj, = 0.9,7 = 1.45, Do = 0.3 X 10 28 cm 2 s _1 (three dots - 
dash). The values of 7 and Do have to be altered in order to ensure that the fits to the radio emission from Fig. 5 
are identical for the different values of fi. The source distribution is parameterized as in Fig. 5. 



region depends therefore sensitively on the assumed inclina- 
tion angle. However, only for a severe underestimate of the 
true inclination angle of NGC 2146 (i.e. if i > 80°) this con- 
clusion cannot be maintained. In this case the radio emission 
along the minor axis would be almost entirely from the halo. 

6.2 CR propagation in the central disk and halo 

The diffusion model is a good fit to the radio data for 
NGC 2146 along the major axis of the galaxy and along the 
minor axis for distances less than 10 arcsec (2.6 kpc) from 
the major axis. This suggests that diffusion is the dominant 
mode of propagation in this region which corresponds to the 
inner, most actively star forming disc of NGC 2146. 

At distances greater than 10 arcsec from the centre of 
the galaxy along the minor axis the model is no longer 
an adequate fit to the data (Figs. 5 and 6). The shape of 
the 1.5 GHz emissivity is very fiat suggesting that trans- 
port of CR electrons is more efficient than would be ex- 
pected from diffusion alone. Furthermore, the spectral in- 
dex flattens which indicates that (i) the energy losses have 
decreased rapidly, or (ii) the diffusion coefficient increases, 
or (iii) a much faster process than diffusion is responsible 
for the transport of the CR electrons. Upon examination of 
the radio emission in the region where the diffusion model 
begins to fail we note that the radio images (Figs. 1 and 3) 
show filamentary structure. The most likely reason for this 
behaviour is that away from the centre of the galaxy a large 
fraction of the radio emission is from a halo which is seen in 
projection. 

The propagation of CR electrons in the halo might be 
more complicated than we have considered for the propa- 
gation within the disc. Evidence for this comes from obser- 
vations of the edge-on galaxies NGC 891 and NGC 4631 
where extended radio continuum emission in the halo with 
a rather flat spectral index has been observed (Hummel 
1991). Breitschwerdt (1994) has explained this flat spectral 
index in a galactic wind model including diffusive and con- 



vective transport of CR's. Siemieniec and Cesarsky (1991) 
modelled the spectral index in the halo of NGC 891 by a 
diffusion model with an outwardly increasing diffusion coef- 
ficient. Both galaxies show that in the halo the propagation 
of CR electrons is faster than that predicted by steady-state 
diffusion alone. 

Due to the strong star formation activity in NGC 2146, 
an outflow from the disc triggered by correlated SN explo- 
sions may be present, a mechanism which is variously de- 
scribed as a "galactic fountain" (Shapiro & Field 1976) or 
"chimney" model (Ikeuchi 1988, Norman & Ikeuchi 1989). 
Such an outflow would transport CR's and magnetic field 
into the halo and could therefore account for a higher ra- 
dio emission in the halo than expected from a pure diffusion 
model. The analysis of X-ray and optical data of NGC 2146 
(Armus et al. 1995) has shown that there are indeed indica- 
tions of the existence of such a starburst-driven superwind. 

6.3 Star formation in NGC 2146 

Hutchings et al. (1990) proposed that NGC 2146 is in the 
late stage of a merger. The merger must have begun > 10 9 
years ago as the rotation curve does not appear to be dis- 
rupted in the outer regions of the galaxy (Young et al. 1988a, 
show the inner regions of the galaxy to have non-circular mo- 
tions of the order of 100 km s~ ). In this scenario material 
has collapsed into the centre of the galaxy and a starburst 
has commenced in the nucleus. 

The data and model presented in this paper further 
emphasize that active star formation is confined only to the 
centre of the galaxy. The radio images show that the star 
formation activity is located in a peculiar spiral S shape 
around the peak of the radio emission. The merger may have 
triggered the formation of a bar which feeds the nuclear star- 
burst with material and is responsible for the unusual struc- 
ture of the radio emission. Star formation is proceeding all 
along the central bar, and the nucleus is not prominent in 
the high-resolution radio images. Young et al. (1988a) sug- 
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gest that if the galaxy is in the late stage of a merger the 
nucleus of the merging companion could be obscured by the 
dust lane. The images of the radio emission presented here 
argue against this as we do not find two spatially distinct 
peaks of radio emission which could be identified as the nu- 
clei of the merging systems. Instead, the star formation is 
distributed throughout the bar-like central region. The un- 
resolved sources seen on our highest resolution images have 
sizes of less than 100 pc and may represented isolated pock- 
ets of star formation or could be individual high-luminosity 
supernovae as seen in M82 (Muxlow et al. 1994). The pres- 
ence of dynamical features such as the arm in the Ha emis- 
sion (Young, Kleinmann & Allen 1988b) and the bar in the 
radio emission suggest that this starburst is dynamically 
driven, i.e. due to a merger or an interaction. 
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